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From the point of view of time-dependent Ginzburg Landau (TDGL) theory, a Josephson-like rela- 
tion is derived for an Abrikosov vortex lattice accelerated and deformed by applied fields. Beginning 
with a review of the Josephson Relation derived from the two ingredients of a lattice-kinematics 
assumption in TDGL theory and gauge invariance, we extend the construction to accommodate a 
time-dependent applied magnetic field, a floating-kernel formulation of normal current, and finally 
lattice deformation due to the electric field and inertial effects of vortex-lattice motion. The result- 
ing Josephson-like relation, which we call an Extended Josephson Relation, applies to a much wider 
set of experimental conditions than the original Josephson Relation, and is explicitly compatible 
with the considerations of TDGL theory. 



I. INTRODUCTION 



The well-known Josephson Relation 1 (JR) relates ap- 
plied electric and magnetic fields to motion of the 
Abrikosov vortex lattice in a thin superconducting sam- 
ple. Taking a simple form, 



(E) 
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the JR is easily understood intuitively via arguments in- 
volving the Faraday Law, and motion of magnetic flux, 
and thought to represent an approximation of general 
applicability, therefore being widely employed by exper- 
imentalists. It is thus time-tested and, though not of a 
fundamental theoretical nature, is an essential and use- 
ful tool. Here, E and B will denote the local electric and 
magnetic fields, and (■ • • ) denotes an averaging in space. 

Relation ([T]) has been known for decades, and does also 
stand up to more modern considerations in the context 
of time-dependent Ginzburg-Landau (TDGL) theory. In- 
deed, the JR is derived by Kopnin^, using some simple 
assumptions about the rigidity of the vortex-lattice mo- 
tion. In this analysis, the premise is that the vortex lat- 
tice already represents some solution to TDGL theory, 
and its motion must remain compatible with that the- 
ory; therefore there is an implicit assumption that the 
vortex lattice is triangular. 

More recently, a derivation along the same lines but ac- 
commodating inertial effects of the condensate has been 
exhibited 3 , which reproduces from the point of view of 
TDGL theory an Inertial Josephson Relation (IJR), and 
which had been obtained much earlier in a hydrodynamic 
context^. The IJR contains a term reflecting the iner- 
tia of the condensate itself, and is thus an extension of 
the JR, applicable at high frequencies. We shall review 
this result below in section IIII1 obtaining the IJR using 
the requirement of gauge invariance, combined with sim- 
ple ingredients of rigid lattice motion and the current 
obtained from TDGL theory. 

In the present work we shall discuss ways in which 
this analysis can naturally be taken further, producing 
what we shall call an Extended Josephson Relation. The 
Inertial Josephson relation mentioned above represents 



one extension of the Josephson relation to accommodate 
condensate inertia. By Extended Josephson Relation, 
we shall mean a relation further extended to include the 
effects discussed below. 

The starting point is chosen to be the Floating- 
Kernel^ (FK) version of the TDGL equation, which 
contains the normal current; so named because it can be 
understood by shifting to the reference frame 'floating' 
with this normal current. 

Not only can a normal-current correction be included, 
but it is not necessary to work with the assumption of 
a rigid Abrikosov vortex lattice. Of course, should the 
(average) magnetic field change with time, so must the 
density of unit-flux vortices. 

Finally, since the electric field and also the inertia of 
the vortex lattice itself provide anisotropic stimuli in the 
plane of the vortex lattice, one ought consider deforma- 
tions of the lattice which include not only the density 
fluctuations mentioned above, but a group of geomet- 
rical deformations which represent a response to these 
anisotropics. 

An Extended Josephson Relation is calculated in sec- 
tion [V] which accommodates all these novel ingredients. 
The result can be considered an extension of the Joseph- 
son Relation (and Inertial Josephson Relation) which can 
describe density fluctuations of the vortex lattice which 
allow for a time- varying magnetic field, and anisotropic 
deformations of the vortex lattice from the electric field 
and vortex acceleration. 



II. NOTATION, CONVENTIONS AND BASICS 

We consider the usual configuration of a planar sample 
in the x-y plane, with E _L z and vortices produced by 
B || z. We will make use of the usual definitions E — 
-\d t A - and B = V x A. 

We write a gauge transformation as 
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where 9 is a function of space and time. Our calculation 
will be gauge- invariant; we will make no gauge choice. 
This allows us to use the principle of gauge covariance. 



III. IJR FROM TDGL 

Within this section, we carefully extract the Inertial 
Josephson Relation from the TDGL equation. We follow 
closely the derivation of the Josephson Relation given by 
Kopnin^ or the IJR given in [3| ; we pay special attention 
to gauge invariance, in anticipation of extending the IJR 
later. 

The TDGL equation is 

D 2 ip + cnp + (3ip\ip\ 2 
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where the gauge-covariant derivative is 



Di/> := (V - tt-AM. 

Tic 



{d t + l -e^+-d t \^ (5) 
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3 S resulting from ([5]) is given by variation of A; 

J, = — ^fD-DW. (7) 



2m* v ' 

Now we must formulate the assumption of rigid motion 
of the vortex lattice. We define by ipo the configuration 
at time zero, ipo(r) := i/j(r,0). With the displacement of 
the vortex lattice given by rz,(t), we thus require 

|V>(r,t)| 2 = hMr-r L (i))| 2 (8) 

so that for some pure-gauge function uj we have 

iP(r,t) = e- iu ^Mr-r L (t)). (9) 

This equation must be gauge- invariant; transforming 
both sides using 9{v,t) (see ((2]) - (j4}), we see that w 
must have the gauge transformation law 

w(r, t) -> w(r, t) + e - (0(r - r L (i), 0) - 9(v, t)) . (10) 

Whatever physical assumption we specify, in the present 
case equation ©, must be independent of any gauge 
choice. Although we do not at this stage know the form 
of the function w, we do know how it must transform, 
and we shall make use of the transformation law (|10p in 
what follows. 

Writing v L = c^r^, we use (|9|) to calculate 

$[$(D-D)V>] = -Vi-V[^(D-D)^] (11) 
+ 2i|^| 2 (v L -V + a t )(Vw + ^-A). 



Equation (IT0|) implies that (vj, • V + d t )oJ gauge- 
transforms as 



(v£-V + ft)w-»-(vL-V + ft)w-— (v L -V + ft)fl, (12) 
which determines 



(vi • V + fltjw = ^(0 - i Vi • A) 



Using 



e*|V>| 2 



(13) 



(14) 



we then obtain 



9 t v s = -v L .Vv s + ^ 7 [iv(v L A)-iv L -VA+E]. (15) 

m* L c c 

Finally, use the identity V(v • V) - v ■ VV = v x V x V 
for constant v to write 

d t v s = -Vi • Vv s + -^-v L x B + ^-E. (16) 
m c m 

Note that this tells us the field v s does not move with 
the lattice; otherwise we would have (<9 t + Vi ■ V)v s = 0. 

We may now take the unit-cell average. Since v s is 
periodic, the first term does not contribute and 



9 t (v s ) = ^ 7 (E) + ^ r v i x(i?). 

m* m*c 



(17) 



This is the IJR as presented in Q. It can be thought of 
as the consequence of describing a rigidly moving vortex 
lattice using TDGL theory, and includes the inertial term 
dt(v s ), absent from the original Josephson Relation. 



IV. FLOATING-KERNEL TDGL 

Our starting point will not be the TDGL equation as 
presented in ([S]), but a version of TDGL supplemented 
by a floating kernel (FK) term^i 



( _ jftv -A J„) V + cal> + /3?/#| 2 

2m* c en 



(18) 



The supercurrent is defined though variation of A, 

J S = ^(|(D-D)-— Jn)V- (19) 
m* 2 en 

In the spirit of a two-fluid model of superconductivity, 
we write 



J s =e*H 2 (v s -v„) 
and write the total current as a sum 

J := J s + J„ 

= e*n s v s + en n v n 



(20) 

(21) 
(22) 
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where v„ := J„/en, J„ = cr„E and n = n n + 2n s . 

In the following section, we do not attempt to solve 
TDGL theory to determine the density or deformation 
dynamics of the vortex lattice, just as in the previous sec- 
tion there was no attempt to solve TDGL to determine 
rigid motion of the vortex lattice. Instead, we allow that 
it may happen and study the consequences. In fact, rigid 
lattice motion may have been too strict an assumption 
even in the case of no normal current; there is no a pri- 
ori reason why the magnetic field must always remain 
constant, and the lattice rigid. 



V. EXTENDED JOSEPHSON RELATION 

In this section we allow much greater freedom for the 
moving vortex lattice; in particular, we allow the magni- 
tude of B to change with time, so that the vortex lattice 
density is time-dependent. Further, we allow the shape 
of the vortex lattice to undergo a global time-dependent 
deformation, in response to anisotropic stimuli. 

We shall generalise the rigid-motion requirement © to 

^(r, t) = e- ffl (''^(A(i)E(i)(r - r L (t)), 0). (23) 

Here A(i) is a dynamic scaling factor for the vortex lat- 
tice. We can expect such a simple dynamic scaling to be 
valid for a B field which does not vary excessively or too 
quickly with time. We require r^(0) = and A(0) = 1, 
and define Bq as the magnitude of the average magnetic 
field when the vortex lattice density is A = 1, therefore 
B(t)^(\B(t)\)=B Q X 2 (i). 

£(£) is a two-dimensional dynamic deformation ma- 
trix for the vortex lattice. Since we have accommodated 
scaling with A, we shall require that £ be an element of 
SX(2, K), and we will write 



and 



E(t) 



= e s(t \ 



(24) 



We shall leave discussion of the matrix S to the following 
section. 

Beginning with the vortex lattice (f2"3"f and the TDGL 
equation with floating kernel (|18p . we proceed as in the 
previous section. Using equation (I23|) we calculate 



where 



T 



d t ip — Y • VV> + iO^ 



Y • Vui — dtu>, 



(25) 



(26) 
(27) 



so that 



d t T>4) = (T • V + ityDip + V(Y • V + iSl)i/> 
ie* 



he 



(Y ■ V — d t )Aijj 



(28) 



. V(D - D)V> rv „ . </>(D - T>U 

d t , , lo = Y-V + VY-- 



2ie* 



he 
- 2ivn. 



IV'I 2 

[dt - Y • V - VY-]A 



(29) 



This allows us to use ([20|) and evaluate 

3*(v s - v„) = (Y • V + VY-)(v s - v„) 
■ [V(Y • A) — Y x B] 



m c 
e 



m 
2 



[E + + — Vfi] 



(30) 



where this time we have used the identity 

V(v • V) - v VV = v x VxV + Vv-V. (31) 

Now, although Y is gauge-invariant, it may be seen from 
(|27p or (|2"5"|) that il must transform as 



fl^n [Y • V0 - d t t 

h 



which determines 



(32) 



(33) 



si = ~[-r-A + (i>]. 

n c 

Substituting ft and Y into (|3U)l . we find 



d A 

d t {v s - v„) = [~7-(r - r L ) • V - v L • V] (v s - v„) 

+ [(r - r L )d t S T V - ^ + d t S T ] (v s - v„) 



+ — i - rd t )E ~ — -r r - rL ^ 

m* me A 

+ d t S(r - r L ) — Vi] x B 



(34) 



where r = m*a n /2e 2 n. Finally we take the spatial aver- 
age and find 

9 t (v s - v„) = ^-(1 - rft)(E) + ^v L x (B) 

+ [^-^]((v s -v n ) + ^„(E)) 

+ ^H|fr L +d t Sr L ] x(B) (35) 
m*c L 2i> J 

where we have substituted for A. We have used that v s is 
periodic and v„ is uniform, and for the sake of simplicity 
we have taken the origin r = to be at the centroid of 
the sample. Equation (|35|) is our Extended Josephson 
Relation (EJR), the main result of this paper, and we 
pause to comment. The first line of Q35p is simply the 
Inertial Josephson Relation of (|17l) . with the addition of 



4 



the r term accounting for the floating-kernel normal cur- 
rent contribution. The remaining corrections consist of 
parts which depend on the time derivative of the mag- 
netic field, and parts which depend on geometrical lattice 
deformation; when deformation is absent, S — 0. 

In the following section we complete the expression of 
the EJR (|35p by showing the explicit parametrisation of 
the deformation matrix S, 



VI. LATTICE DEFORMATION 

A two-dimensional space may be expanded by a scale 
factor e CT in the x direction, and e~ a in the y direction 
by the matrix 



£(0,<r) 



e a 
e~ c 



(36) 



which has unit determinant. Rotating so that the ex- 
pansion is along a line at angle $ and the contraction 
perpendicular, 



£(0,<7) 



fl(-tf)£(0, o-)-R(tf) 

' 1 - 2sin 2 tf 2sinz9cosz9 
2 sin -d cos i9 1 — 2 cos 2 i? 



cxp < a 



exp 5 P 



(37) 



Here we have parametrised the deformation by fi := 
(crcos$,crsin$). This set of deformations is not closed; 
composition can produce rotations. When combined with 
rotations all elements of SX(2,K) can be constructed. 
Some identities are E^cr) -1 = £(i? + 7r/2,er) = 
£(i?,-cr), £(0,0) = 1 and E(i?,02)£(i9,oi) = £(0,oi + 

If we define convenient matrices which depend on the 
deformation parameter /i, 



Ma? My 

My Ms 



and 



Ma; My 
My M^ 



(38) 



we may write Su, = w„n„/ /x. For a given time-dependent 
deformation parameter we may calculate 



"'u r 

2n„/ 

M 



M 4 M 

M 2 



(39) 



Now, let us recall the 5 matrix introduced in section 
IVl equation (|24l) . Given some deformations parametrised 



by M; 



we would set S 



In fact, 



it is possible immediately to write down several plausible 
examples of physical sources of deformation. 

Let us anticipate that the application of an electric 
field (in the plane) , or a time-derivative of this field may 
tend to deform the vortex lattice. We can accommodate 
this with a deformation parameter 



(40) 



Another possibility is that the global inertia or accelera- 
tion of the vortex lattice would coincide with a deforma- 
tion; in this case, 



M : = vqv l + vid t VL- 



(41) 



Obviously, there is an approximation involved here in 
that our consideration is limited to a global deformation. 

Adding the above two deformation contributions to- 
gether, 



S(t) 



(42) 



dtS(t) may be calculated for each of the two terms by 
using (]3gj); 



d t S(t) 



M 

TO,, 



2n, 



2n,v 



MM 

M 2 
v ■ v' 



(43) 



Upon substitution of dtS and dtS T the Extended Joseph- 
son Relation (|35p will contain the parameters Eq, E±, vq 
and v\. In principle, these parameters could be fit to 
experimental data, characterising response beyond the 
usual JR or IJR. 



VII. CONCLUSIONS 

As we mentioned in the introduction, the Josephson 
Relation, though simple in form, is often used by exper- 
imentalists to understand the motion of vortices in the 
presence of an applied field. It has been shown in the 
literature how to recover the Josephson Relation^ and 
a form valid at higher frequencies, the Inertial Joseph- 
son Relation^ in the context of TDGL theory. We have 
taken this technique further, using the assumption of a 
vortex lattice solution to TDGL theory, to extend the 
Josephson Relation to a form covering a floating-kernel 
formulation of TDGL with a normal current in the spirit 
of a two-fluid model. Additionally, we have shown how 
to include effects of time- varying magnetic flux, and also 
global vortex-lattice deformations; our considerations are 
expected to be valid for small changes in magnetic field 
and small lattice deformations. This is not due to any 
approximation in the calculation itself, rather it is due to 
the implicit assumption of the characteristics of the vor- 
tex lattice dynamics, equation (|23p , which nevertheless 
represents a far weaker assumption than that of perfectly 
rigid global motion ©. 
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